Uniform spin chain physics arising from NCN bridges in CuNCN: 
surprises on the way from copper oxides to their nitride analogs 
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We report on the unexpected uniform spin chain physics in CuNCN, the insulating nitride analog 
of copper oxides. Based on full-potential band structure calculations, we derive the relevant micro- 
scopic parameters, estimate individual exchange couplings, and establish a realistic spin model of this 
compound. The structure of CuNCN contains chains of edge-sharing CuN4 squares. As a surprise, 
in contrast to analogous [CUO2] chains in "edge-sharing" cuprates, the leading magnetic interactions 
J ~ 2500 K run perpendicular to the structural [CUN2] chains via bridging NCN groups. The result- 
ing spin model of a uniform chain is in agreement with the experimentally observed temperature- 
independent magnetic susceptibility below 300 K. The nearest-neighbor and next-nearest-neighbor 
interactions along the structural [CUN2] chains are Ji ~ —500 K and J2 ~ 100 K, respectively. De- 
spite the frustrating nature of Ji and J2 , we assign the susceptibility anomaly at 70 K to long-range 
magnetic ordering, which is likely coUinear with antiparallel and parallel arrangement of spins along 
the c and a directions, respectively. The pronounced one-dimensionality of the spin system should 
lead to a reduction in the ordered moment and to a suppression of the transition anomaly in the 
specific heat, thus impeding the experimental observation of the long-range ordering. Our results 
suggest CuNCN as a promising material for ballistic heat transport within spin chains, while the 
sizable bandwidth W ~ 3 eV may lead to a metal-insulator transition and other exotic properties 
under high pressure. 
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I. INTRODUCTION 

One-dimensional (ID) electronic systems are in the fo- 
cus of experimental and theoretical research due to the 
exotic properties that emerge in low-dimensional models 
and can be observed in real systemsi^ Coppcr(II) ox- 
ides with chain-like structures are one of the best play- 
grounds for studying ID physics of localized electrons, 
because the strong Coulomb repulsion in Cu d shell and 
the half-filling regime lead to insulating behavior, while 
the rich crystal chemistry of copper oxides allows to vary 
relevant microscopic parameters. Copper oxides can be 
properly described within simple or extended Hubbard 
models."^ Subsequently, the low-energy properties of un- 
doped (Cu+^-containing) systems are easily mapped onto 
the spin-only Heisenberg model4ii^ 

The two main scenarios of Cu-based ID systems are 
the so-called " corner-sharing" and " edge-sharing" chains 
that correspond to corner or edge connections between 
adjacent CUO4 plaquettes within the chain. In "corner- 
sharing" chains, the Cu-O-Cu angle is usually close to 
180° and allows for the strong superexchange between 
neighboring Cu atoms. This type of the chain struc- 
ture is found in Sr2Cu03, the archetypal material for 
the Heisenberg model of uniform spin-^ chaiu)^ This 
model allows for ballistic heat transport mediated by 
spin excitations.®'^ Indeed, Sr2Cu03 and other spin-chain 
materials reveal surprisingly high thermal conductivities 
along the direction of spin chainsj^^ 

The case of the "edge-sharing" chains is quite differ- 
ent. Copper atoms are located closer to each other, 
and two competing magnetic interactions emerge. The 
nearest-neighbor (NN) coupling Ji is usually ferromag- 



netic (FM), because the Cu-O-Cu angle is close to 90°. 
In contrast, the next-nearest-neighbor (NNN) coupling 
J2 is antiferromagnetic (AFM) and corresponds to Cu- 
0-0-Cu superexchange. The absolute values of Ji and 
J2 are usually comparable. This leads to strong mag- 
netic frustration. For J2/J1 < —0.25, the energy of the 
system can be minimized in the spiral state, where ad- 
jacent spins are turned for a constant angle with respect 
to each other. Although the Mermin- Wagner theorem 
states the lack of the long-range ordering in a purely ID 
system above zero temperature, non-negligible interchain 
couplings stabilize the finite-temperature spiral order in 
real materials. This type of the ordering was recently ob- 
served in a number of "edge-sharing" compounds, such 
as LiCu202,iiii^'i^ LiCuV04,^'''^^ and Li2CuZr04 J^ii^ 

A spiral magnetic ground state discloses further un- 
usual phenomena. Below the ordering temperature, 
LiCu202 and LiCuV04 show unconventional ferro- 
electricity which is strongly coupled to the magnetic 
fieldi ^^i^^ However, the origin of this behavior remains 
controversial. The ferroelectricity can be explained ei- 
ther as a purely electronic effect (the tendency to re- 
lease the frustration)^ or as a result of Li/Cu anti- 
site disorderi2ii2^ To get further insight into the fer- 
roelectricity of Cu-based spin-chain materials, one has 
to study other compounds showing the spiral ground 
state, e.g., with different ligands. The size, charge, and 
chemical nature of the ligand control the relevant micro- 
scopic parameters: the hoppings and the magnitude of 
electronic correlations (the Coulomb repulsion in Cu d 
shell, which is partially screened by the ligand orbitals). 
The search for new materials led to the recent study of 
CuCl2 with chains of edge-sharing CUCI4 squaresi^i^ 



2 



Low-temperature neutron diffraction evidences its spiral 
magnetic structure^i and calls for the further study of 
the possible low-temperature ferroelectricity in this com- 
pound. Another suitable ligand is nitrogen. However, ni- 
tride compounds of Cu^^ are scarce, because the nitride 
ligand usually stabilizes low oxidation states of transi- 
tion metals - e.g., -1-1 for Cu in the simple copper nitride 
CU3N. 

Recently, the nitride environment of Cu^'^ was 
achieved in copper cyanodiimide CuNCN.^^ The crystal 
structure of this compound reveals chains of edge-sharing 
CuN4 plaquettes arranged along the a axis. The nearly 
linear NCN groups couple the chains in the ac plane, 
while the resulting layers are stacked along the b axis 
(Fig. [T|). CuNCN has black color. Resistivity measure- 
ments on polycrystalline samples evidence semiconduct- 
ing behavior with the activation energy of 0.1 eV. Below 
300 K, the magnetic susceptibility of CuNCN is nearly 
temperature-independent. Additionally, a susceptibil- 
ity anomaly is observed at 70 K, whereas no respective 
anomaly is seen in the specific heat. The low-temperature 
neutron diffraction does not manifest long-range mag- 
netic orderingi^ 

Band structure calculations suggested a spin model of 
an anisotropic triangular lattice with leading exchange 
couplings in the ab plane. However, the relevant mi- 
croscopic parameters of the Hubbard model were not 
derived. Additionally, the proposed spin model seems 
to be counter-intuitive from chemical point of view, be- 
cause the strong interaction along the b direction {Jb of 
about 800 K according to Ref. [2^, see Fig. |4] for no- 
tation) corresponds to the rather long Cu-Cu distance 
of 3.43 A and lacks any obvious superexchange pathway 
(see Fig. [T]) . The reference to oxide materials with " edge- 
sharing" chainsiiii^ii^ii^ suggests that the half-filled Cu d 
orbital lies in the CUN4 plane, while four other d orbitals 
are fully occupied and do not take part in the magnetic 
interactions. Then, the interlayer couplings should be 
negligible, whereas strong exchange couplings should run 
along the structural chains and show magnetic frustra- 
tion. 

In the following, we perform accurate full-potential 
band structure calculations that evaluate the relevant mi- 
croscopic parameters and suggest a realistic microscopic 
model for CuNCN. We show that the "on-site" physics of 
this nitride material indeed resembles copper oxides, the 
only difference being a reduction in the on-site Coulomb 
repulsion due to the stronger screening caused by ni- 
trogen as a ligand. However, the "inter-site" physics is 
dramatically changed, because NCN groups mediate un- 
expectedly strong superexchange couplings. Then, the 
magnetic behavior is ID. However, the spin chains do 
not match the structural chains but rather run perpen- 
dicular to them. We carefully compare our results to the 
available experimental data, analyze similarities and dif- 
ferences between CuNCN and copper oxides, and provide 
an outlook for further experimental studies. 




FIG. 1: (Color online) Crystal structure of CuNCN showing 
chains of edge-sharing CUN4 plaquettes. The chains are linked 
into layers via NCN groups. The Cu atoms are located in the 
centers of the green CUN4 plaquettes. 



II. METHODS 

Our microscopic model is based on band structure cal- 
culations performed within local density approximation 
(LDA) of density functional theory. We use the exchange- 
correlation potential by Perdew and Wang^"^ and the full- 
potential band structure code with a basis set of atomic- 
like local orbitals (FPLO8.50-32),28 To check the robust- 
ness of the results, we performed a number of compara- 
tive calculations applying the generalized gradient ap- 
proximation (GGA) with the Perdew-Burke-Ernzerhof 
exchange-correlation potential.— In order to compare our 
results to the previous study (Ref. [1^, we also repeated 
several calculations using the projector augmented wave 
(PAW) method'^*' implemented in the Vienna ab-initio 
simulation package (VASP).— 

For the LDA calculation, the proper partitioning of 
the reciprocal space is achieved by the k mesh with 2041 
points in the irreducible part of the first Brillouin zone. 
The convergence with respect to the k mesh was care- 
fully checked. We used two sets of structural parameters 
derived from the refinement of the room-temperature x- 
ray data3^ and the low-temperature neutron datai^ The 
two sets of atomic coordinates yield essentially similar 
results, hence one should not expect any pronounced 
temperature-dependent changes in the electronic struc- 
ture of CuNCN. 

While LDA provides a simple picture of electronic 
states and yields input parameters for the free-electron 
part of the Hubbard hamiltonian, it misses the strong 
electronic correlations in the d shell and fails to de- 
scribe the correlation-induced band gap in many transi- 
tion metal compounds. To account for the correlation 
effects, one has to include them either on the model 
level (on top of self-consistent band structure calcula- 
tions) or into the self-consistent procedure itself. For 
the model treatment of correlations, we fit relevant LDA 
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bands with a tight-binding (TB) model and include the 
resulting hopping parameters into a one-orbital Hubbard 
model with the effective on-site Coulomb repulsion [/off- 
Then, the Hubbard model can be solved using numeri- 
cal simulations. For the low- lying (spin) excitations and 
for the half-filling regime, we can reduce the hamilto- 
nian to a Heisenberg model and calculate the exchange 
integrals via the simple expression of second-order per- 
turbation theory. Further details of the procedure are 
given in Sec. IIIII 

To include electronic correlations into the self- 
consistent procedure, one usually applies the local spin- 
density approximation (LSDA)-t-C/ approach that treats 
correlations in a mean-field way. LSDA+C/ gives reli- 
able energies for different types of spin order, and these 
energies can be further mapped onto the classical Heisen- 
berg model. This approach has already been utilized in 
Ref. I26. However, the authors of Ref. [2^ performed cal- 
culations for a relatively small supercell which allowed to 
access few exchange couplings only. In the following, we 
will show that further, long-range couplings should also 
be included in the model. To evaluate them, we use two 
32-atom supercells: 2a x b x c and 4a x i(a + b) x c, 
where a, b, and c are the translation vectors of the crys- 
tallographic C-centered unit cell of CuNCN. The k mesh 
included 256 and 192 points, respectively. The result- 
ing total energies allow to estimate all the couplings ev- 
idenced by the LDA-based analysis. The use of the two 
supercells has an additional advantage, because short- 
range couplings are independently estimated in each of 
the supercells. Then, the comparison of the resulting 
numbers suggests a natural error bar for the computa- 
tional results. In contrast to other copper compounds 
(e.g., Ref. |33 ). the two supercells yield slightly different 
estimates of the exchange couplings. This issue will be 
further discussed in Sec. IIIII 

The electronic correlations are parametrized by the ef- 
fective on-site Coulomb repulsion potential UeS of the 
Hubbard model or by the Ud (repulsion) and Jd (ex- 
change) parameters of LSDA-I-C/. To get an idea about 
the values of these parameters, we use the constrained 
LDA approach^^ implemented in the TB-LMTO-ASA 
(tight-binding linearized muffin-tin orbitals in atomic 
spheres approximation) code.'^* Under the constraint of 
zero hoppings between correlated and uncorrelated or- 
bitals, we calculate the energy of the correlated (Sd) 
orbital for several fixed occupancies. This yields the 
estimates of t/^MTO ^ 5.6 eV and jLmto ^ ^^.0 eV 
for CuNCN. These numbers can be further used for 
LSDA-hf/ calculations within LMTO-ASA. However, the 
transfer to other codes or models requires a more care- 
ful consideration of the meaning regarding [/LMTO ^^^^ 

jLMTO 

The constrained LDA procedure yields the Ud and Jd 
parameters which are relevant for the specific (muffin- 
tin) orbitals of the LMTO code. In other codes, these 
parameters are applied to different d functions, hence 
their values should also be different. In the Hubbard 



model, one treats correlated bands rather than atomic 
orbitals. In the case of CuNCN, nearly half of the corre- 
lated bands originate from uncorrelated nitrogen states 
(see Sec. IIIII) . hence a strong screening is expected, and 
the UcS value should be reduced compared to Ud and 
Jd- Thus, the [/LMTO ^^^^ jLMTO numbers can only 
be used in a comparative way and can not be trans- 
ferred to LSDA-|-[/ within FFLO or to the Hubbard 
model treatment. For the proper comparison, we per- 
form a constrained LDA calculation for a similar copper 
oxide material LiaCuOa.^^ We find [/LMTO ^ 9 4 and 
jLMTO ^ i Q gv. The jLMTO ^j^j^g ^j^gg depend 
on the ligand, while [/^mto ^^j. Li2Cu02 is strongly en- 
hanced compared to [/LMTO ^ 5.6 eV for CuNCN. This 
analysis suggests that the on-site Coulomb repulsion is 
ligand-dependent. Nitrogen p orbitals show larger over- 
lap with Cu d orbitals, thus leading to improved screening 
of the on-site Coulomb repulsion. Thus, nitride materi- 
als will generally require lower U values compared to the 
respective oxides. 

After establishing this trend, we can turn to our ex- 
perience of band structure calculations for copper ox- 
ides. For copper oxides treated within FPLO, the 
Ud values in the range from 6.0 to 8.0 eV are now 
well established^i26i37^ ^he comparison of [/LMTO ^j. 

CuNCN and LiaCuOa (6.6 eV vs. 9.2 eV) suggests lower 
L/d of 5 - 6 eV for CuNCN in FPLO. For Jd, we fix 
the value of 1.0 eV. In the case of the one-orbital Hub- 
bard model, the Ucs values of 3.5 — 4.5 eV have been ob- 
tained for copper oxides by fitting model simulations to 
different experimental results (exchange couplings, opti- 
cal conductivity, electron energy loss spectra, etc.)^'^'^! 
For the nitrides, it is reasonable to take a somewhat lower 
UcS, and we use UcS = 3.5 eV in the present analysis. 
Further experimental studies of CuNCN should provide 
a more accurate estimate of this parameter. 



III. RESULTS 

A. LDA-based analysis 

We start with the analysis of LDA band structure of 
CuNCN. The density of states plot is shown in Fig. [51 
Valence bands of CuNCN are formed by Cu 3d and N 
2p orbitals with the minor contribution of carbon states. 
Empty bands above 2 eV correspond to antibonding 
states of the NCN groups. The energy spectrum is metal- 
lic, because LDA does not account for correlation effects 
and underestimates band gaps in transition-metal com- 
pounds. The energy gap can be reproduced in LSDA+[/ 
(see below). 

In CuNCN, the local environment of copper is subject 
to a strong Jahn- Teller distortion. We find four short 
Cu-N distances of 2.00 A within the CUN4 plaquettes 
and two longer Cu-N distances of 2.61 A. Similar to 
oxides, this type of the local environment should lead 
to a pronounced crystal field splitting with the highest- 
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FIG. 2: (Color online) LDA density of states for CuNCN. The 
Fermi level is at zero energy. 
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FIG. 3; (Color online) LDA bands (thin black lines) and 
the fit of the TB model (thick green lines). Dots show 
the contribution of the Cu d^2_y-2 orbital. The notation 
of k points is as follows: r(0,0,0), X(x,0,0), S(f,0.5, 0), 
F(0, 0.5,0), Z(0, 0,0.5), 5(0,0.5,0.5), i?(f , 0.5, 0.5), and 
T{x, 0, 0.5), where x = 0.25 + a^/4b^ ~ 0.3087, and the coor- 
dinates are given along kx,ky, and in units of the respective 
reciprocal lattice parameters 4n /a, in /b, and 2n/c. 

lying orbital having — character (x axis runs along 
one of the short Cu-N bonds, while the z axis is per- 
pendicular to the CuN4 plaquette)^i^i2a The respec- 
tive orbital character is easily recognized in the LDA 
bands near the Fermi level (Fig. [3]). We find two dj;2_y2 
bands corresponding to two Cu atoms in the primitive 
cell of CuNCN. The hybridization to other bands is weak. 
Therefore, the fit with an effective one-orbital TB model 
is straightforward. To extract hopping parameters, we 
calculate overlap integrals of maximally localized Wan- 
nier functions centered on Cu sites i^i^ The resulting 
hoppings ti are listed in Table HI 

The leading hopping runs along the c direction between 
the structural chains. This hopping t ~ —0.384 eV is very 



TABLE I: Leading hoppings of the TB model U (in eV) and 
the resulting AFM contributions to the exchange integrals 
Jf"™ = U^/U,B (in K), UcB = 3.5 eV. 



t 


ti 


t2 




-0.384 


0.031 


0.058 


-0.073 


jAPM 


jAPM 


tAPM 


tAFM 


1960 


13 


45 


71 



large and leads to the sizable bandwidth W ~3 eV. Still, 
the largest t is one order of magnitude smaller than Uefi, 
hence the perturbation treatment of the corresponding 
Hubbard model should be reasonable. For the case of the 
half- filling and for the low- lying excitations, the effective 
one-orbital Hubbard model is reduced to a Heisenberg 
model with AFM exchange J^™ = itf/UcS- Using this 
simple expression and Ucs = 3.5 eV, we find the AFM 
contributions to the exchange couplings, which are listed 
in Table [J The strongest AFM coupling J ~ 1970 K is 
mediated by NCN groups (see Fig.[4|). The NNN coupling 
within the structural chains (J2) and the long-range cou- 
pling in the ac plane (Jaci) are much weaker and amount 
to 50 — 70 K only. Further hoppings are below 0.03 eV, 
i.e., the respective J^^^ do not exceed 10 K. In partic- 
ular, the AFM interaction along the b direction is negli- 
gible {tt = 0.009 eV, i.e., J^™ ~ 1 K). 

It is also instructive to consider the shape of the Wan- 
nier functions. In Fig. [51 we plot Wannier functions for 
CuNCN and for the typical "edge-sharing" oxide mate- 
rial Li2Cu02.'^^ Each of the Wannier functions is com- 
posed of atomic Cu Sd,j.2_y2 orbital and the 2p orbitals 
of the neighboring nitrogen/oxygen atoms. However, 
the Wannier function for CuNCN has additional con- 
tributions from second-neighbor nitrogen atoms in the 
NCN groups. The reason for this difference is the strong 
TT-bonding within the NCN group. In oxides, Cu or- 
bitals usually overlap with individual (atomic) oxygen 
orbitals. In CuNCN, the Cu dx2_y2 orbital overlaps with 
the highest occupied molecular orbitals (HOMO's) of the 
NCN groups. For the idealized linear NCN group, the 
HOMO is composed of the two nitrogen orbitals and 
has TTg symmetry4i Then, this HOMO enters the Cu- 
based Wannier function and brings significant contribu- 
tion from the second-neighbor nitrogen atoms. The spa- 
tial extension of the Wannier functions illustrates the 
large hopping along the c direction and explains the lead- 
ing magnetic interactions via the NCN groups. 

The one-orbital Hubbard model yields AFM couplings 
only. To estimate FM contributions, one has to extend 
the model and to include nitrogen p states. The largest 
FM contribution is expected for Ji due to the 90° Cu- 
N~Cu superexchange. For this scenario of edge-sharing 
copper plaquettes, Mazurenko et al^ have recently pro- 
posed a multi-orbital expression for the total exchange. 
In a nutshell, their result for the exchange coupling be- 
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FIG. 4: (Color online) Magnetic interactions in the ac (left) 
and ah (right) planes of the CuNCN structure. Bottom figures 
sketch the respective projections of the spin lattice. 

tween sites i and j can be written as follows 

J^J = 77^ " 2/3VpiVz, (1) 

f^eff 

where f3 shows the contribution of each ligand to the cor- 
related bands/Wannier functions, Jp is Hund's coupling 
at the ligand site, and Ni is the number of ligands where 
the Wannier functions of adjacent copper sites overlap. 
Compared to the result of Ref. H^, we made two adjust- 
ments of the Eq. 11]). First, we consider Jij as the full 
exchange coupling for the i — j bond, hence both terms 
in Eq. ([T]) are multiplied by a factor of 2. Second, we put 
[/eff in the denominator of the first term (see footnote 
for further details). "''^ These adjustments provide trans- 
parent physical meaning of the Eq. ([T|). The first term 
is the usual AFM coupling resulting from the one-orbital 
Hubbard model. The second term is the energy gain for 
the FM configuration due to the Hund's coupling on the 
ligand site. 

To find /3, we consider the contributions of nitrogen 
atomic orbitals to the maximally localized Wannier func- 
tions centered on Cu sites. Each nitrogen atom in the 
CuN4 square shows ~ 0.088. Assuming the reason- 
able value of Jp = 1.5 eV,i2 we find jf^ ~ -550 K. 
Thus, there is a strong FM coupling along the structural 
chains. In the following, we will confirm this result using 
LSDA-t-t/ calculations. 



B. LSDA+[/ 

According to Sec. [Hi Ud = 5 — 6 eV is a reason- 
able range for the LSDA-I-C/ Coulomb repulsion param- 
eter in CuNCN. Such Ud values are comparable to the 
LDA bandwidth W ^ 3 eV, hence one can consider 
CuNCN as a material on the border between the regimes 



of strong and moderate correlations. To perform cor- 
rect modeling, one can approach the problem from both 
sides: either perform LSDA-I-C/ calculations assuming 
that the correlations are sufficiently strong or solve the 
multi-orbital Hubbard model using dynamical mean-field 
theory (DMFT) approach. In the following, we con- 
sider the first method and leave a DMFT treatment for 
future investigations. The reasons for this choice are 
twofold. First, it is known that LSDA-I-C/ removes the 
quasiparticle band and severely fails in describing spec- 
tral functions and metal-insulator transitions. However, 
LSDA4-C/ should still work for the phenomena related to 
localized electrons - e.g., for magnetic interactions in the 
insulating material CuNCN. Second, weak exchange cou- 
plings can be evaluated at sufficiently low temperatures 
only, while such temperatures are still hardly reachable 
in the widely-used Monte-Carlo solvers of the DMFT im- 
purity problem. Therefore, we restrict ourselves to the 
LSDA-f C/ analysis. We will show that LSDA-I-C/ yields 
reasonable results and seems to work quite well even ap- 
proaching the regime of moderate correlations. 

The proximity of Ud to the LDA bandwidth W leads 
to small energy gaps (Eg) and even impedes the con- 
vergence of the calculations once the gap approaches 
zero. This problem is strongly dependent on the double- 
counting correction and on the spin configuration. We 
will consider both issues in more detail. The double- 
counting correction is a necessary part of the LSDA+C/ 
approach, because the correlation energy is partially in- 
cluded in LSDA and has to be subtracted, once the ex- 
plicit Hubbard-like term is added to the energy func- 
tional. The two main approaches to the double-counting 
correction scheme are around-the-mean-field (AMF)^ 
and the fuUy-locahzed-hmit (FLL).'^''''^'^ In the AMF 
scheme, one derives the double-counting term assuming 
averaged occupancies for all the correlated orbitals. In 
contrast, the FLL approach uses integer (0 or 1) occu- 
pancies. It is generally believed that the AMF scheme 
should be used for materials with moderate electronic 
correlations, while the FLL double-counting correction is 
more appropriate for strongly correlated (fully localized) 
systems4^ However, the role of the double-counting cor- 
rection has been scarcely studied. Recent calculations 
for model systems established the strong underestimate 
of the magnetic moment on rare-earth atoms in AMF 
due to the lack of the Hund's coupling on the correlated 
site.'*'' This, however, is irrelevant to the present situa- 
tion, because copper has one unpaired electron only, and 
any type of the LSDA-I-C/ functional will generally favor 
the magnetic rather than the non-magnetic state. 

For the AMF double-counting correction, spin config- 
urations with AFM ordering along the c axis show lowest 
energies, largest Eg values, and can be stabilized down 
to C/d = 3 eV. For the alternating FM-AFM ordering 
along the c axis, the convergence is reached at Ud> 5 eV 
only. Finally, the FM ordering along the c axis leads to 
the highest energy and reaches the gapped ground state 
at C/d > 7 eV only. To evaluate J, one has to use at 
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FIG. 5: Maximally localized Wannier functions for CuNCN (left panel) and Li2Cu02 (right panel). The dashed line shows the 
CuN4 (CUO4) plaquette. 



TABLE II: LSDA+t/ (AMF) estimates of the exchange inte- 
grals (in K), Ud is the Coulomb repulsion parameter (in eV). 
The notation of the exchange couplings is shown in Fig. |4] 
Further couplings are below 5 K. 





J 


Ji 


J2 


Jac 


Jb 


5 


2900 


-500 


120 


32 


4 


6 


2340 


-410 


98 


52 


2 


7 


1880 


-330 


77 


70 


2 


9 


1210 


-180 


46 


65 


2 



least two of these configurations. Therefore, the natural 
lower limit for calculating exchange integrals in CuNCN 
is J7d = 5 eV. In the following, we will also consider higher 
Ud values for the sake of comparison. 

In CuNCN, the AMF double-counting correction 
scheme shows superior results compared to FLL, because 
the latter yields smaller energy gaps and fails to con- 
verge the sufficient set of different spin configurations for 
Ud<& eV. For Ud>Q eV, FLL and AMF yield similar 
results with a nearly constant offset of 1.5 — 2.0 eV in 
Ud (to reproduce the AMF results with FLL, a higher 
Ud value should be used). This means that the result- 
ing exchange couplings are robust with respect to the 
LSDA+C/ implementation, although the application of 
FLL requires additional caution. We should note that 
the application of LSDA-I-J7 to moderately correlated sys- 
tems has already been tested by Petukhov et al^ They 
used LSDA-f [/ for calculating the band structure of the 
correlated metal FeAl, found better results for AMF, and 
suggested the interpolating scheme between the AMF 
and FLL double-counting corrections. In the present 
case, this scheme would likely lead to a mere shift of 
Ud: because AMF and FLL yield similar results. 

In Table [Til we summarize the exchange couplings ob- 
tained from AMF LSDA-|-[/ calculations. Some of the 
couplings were estimated in two different supercells (see 
Sec. [n]), and different results were found. For example, 
at [/rf = 5 eV the calculations for the 2a x b x c supercell 
yield Ji = —544 K and Jf, = 4 K, while the calculations 
for the 4ax i(a+b) xc supercell lead to Ji+Jb = —454 K. 



The J I = —500 K value in Table HIl is a rough average of 
these estimates. In a similar way, we find 10 % error bar 
for Jac- Upon increasing Ud, the difference between the 
two supercells is reduced. At iJd = 9 eV, both supercells 
show consistent estimates of all the exchange couplings 
similar to strongly localized systems (e.g., Ref . i32[) . This 
feature may be related to the regime of moderate cor- 
relations in CuNCN. 4^ It is also worth to note that the 
change in the exchange-correlation potential (GGA in- 
stead of LDA for one of the supercells) leads to a similar 
10 % error bar for J. This error bar persists even at high 
Ud. 

The above considerations show that the LSDA-I-J7 esti- 
mates of the exchange couplings include the sizable error 
bar of about 10 %. Nevertheless, they yield a robust 
microscopic scenario which is consistent with the LDA- 
based analysis. The leading coupling J runs along the 
c axis and reaches nearly 3000 K assuming Ud ~ ^ eV. 
The NN coupling within the structural chains is indeed 
FM: Ji ~ -500 K. There is also the NNN AFM cou- 
pling within the structural chains (J2 — 100 K) and the 
AFM coupling in the ac plane J ad = 30 — 50 K. Further 
couplings are very weak in agreement with the results of 
the TB analysis. At higher Ud, the scenario persists, but 
the absolute values of the leading exchange couplings are 
reduced. 

Our model analysis and our LSDA-|-(7 calculations con- 
sistently point to the leading exchange couplings along 
the c direction in CuNCN. The visualization of the Wan- 
nier functions suggests the crucial role of the NCN groups 
in mediating these interactions. We also find the siz- 
able FM coupling within the structural chains and weaker 
AFM couplings. These results will be compared to the 
experimental findings in the next section. In the remain- 
der of this section, we will try to find out why our ID 
model with the leading coupling along the c direction is 
different from the previously proposed 2D model with 
leading couplings in the ah planci^ 

While our study and the previous investigation by 
Liu et al^ use a similar approach to evaluate the ex- 
change couplings (mapping LSDA-|-[/ total energies onto 
the classical Heisenberg model), there are several differ- 
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ences which might be responsible for the difference in the 
resulting scenario. These differences are: i) exchange- 
correlation potential (LDA vs. GGA); ii) basis set (lo- 
cal orbitals in FPLO vs. PAW in VASP); iii) double- 
counting correction schemes (conventional AMF and FLL 
vs. FLL-like Dudarev's approach"*^ in VASP). The first 
issue has been discussed above: LDA and GGA yield sim- 
ilar results within FPLO. To test the two other issues, we 
performed comparative calculations in VASP using both 
the conventional FLL and the Dudarev's schemes for the 
double-counting correction (in VASP, the AMF scheme 
is not available). At C/d = 10 eV and Jd — I eV (this 
should roughly correspond to Ud — 9 eV and Jd = 1 eV 
in FPLO), we arrived to the FPLO results within the 
same 10 % error bar. However, the decrease oi Ud down 
to 7 eV (corresponding to Ud — Jd = d eV in Ref. [1^ 
or to Ud = 6 eV in FPLO) dramatically worsened the 
convergence. This, however, is hard to detect, because 
VASP normally achieves the convergence in total energy 
and does not try to reach good charge convergence. At 
Ud = 7 eV, some of the spin configurations lacked the 
energy gap and failed to converge the charge below 10~^. 
Then, the reliability of the resulting total energies might 
be questionable. Moreover, the metallic regime does not 
allow to map the energies onto the Heisenberg model. 
Thus, we can suggest that for moderately correlated sys- 
tems VASP should be used with caution. In this situa- 
tion, band structure codes implementing the AMF ver- 
sion of LSDA+U are preferable. 



IV. COMPARISON TO THE EXPERIMENT 
AND DISCUSSION 

According to our estimates of the exchange couplings 
(see Tables U and [TT]), the uniform spin chain should 
be a reasonable first approximation to the spin model 
of CuNCN. The magnetic susceptibility of the uniform 
spin-i chain shows a maximum at Tmax = T/J ~ 0.6 
with the reduced susceptibility x* = 0.1469, where 
X* — xJ/{NAg'^tJ''B)^ Na is Avogadro's number, g is 
the (^-factor, and hb is Bohr magneton.^" Assuming 
J ~ 2500 lOi and g = 2, we find T^ax 1500 K which is 
well above the experimentally studied temperature range 
between 1.85 K and 320 K.^^ Below Tmax, the suscep- 
tibility decreases down to the finite value Xo — 0.102 
at T = (at zero temperature, the susceptibility re- 
mains finite due to the strong quantum fluctuations in 
a ID system).— The range of 1.85 — 320 K corresponds 
to < T/J < 0.15, where temperature dependence 
of the susceptibility is weak. Indeed, the experimental 
data show nearly temperature-independent susceptibil- 
ity of about 9 • 10~^ emu/mol.^^ This value is in good 
agreement with xo — 6.1 • 10~^ emu/mol for the uniform 
spin-i chain with J ~ 2500 K. The discrepancy can be 
attributed to additional temperature-independent contri- 
butions (core diamagnetism. Van Vleck paramagnetism) 
and to the possible impurity contribution in Ref. 26. To 



get an accurate experimental estimate of J, one has to 
measure the susceptibility up to Tmax (or, at least, sufR- 
ciently close to Tmax). 

ID spin systems do not undergo long-range ordering 
down to zero temperatures, but non-negligible interchain 
couplings usually cause magnetic ordering in real low- 
dimensional materials. In the case of CuNCN, we find a 
sizable interchain interaction Ji that favors FM ordering 
pattern along the a direction. Such a pattern is further 
stabilized by Jac but destabilized by J2 (see Fig. U]). To 
get an idea about the resulting magnetic structure, it 
is instructive to consider the couplings along the a di- 
rection within the frustrated spin chain model. In this 
model, the ordering is FM at J2/J1 > ac = —0.25, while 
at lower J2 / Ji values the spiral ordering is stabilizedi^ 
In CuNCN, Ji ~ -500 K and J2 ^ 100 K, i.e., J2/J1 ~ 
—0.2 is close to ac- On the other hand, the ac value 
should inevitably be modified by the strong coupling 
J > I J2. This regime has not been studied theoreti- 
cally. Nevertheless, one can suggest that the coupling J 
will drive the system away from the spiral ordering, be- 
cause strong quantum fluctuations in the quasi- ID sys- 
tem will favor the collinear ordering via the " order-from- 
disorder" mechanismi^ The relevant experimental exam- 
ple is Li2Cu02, where the J2/J1 ratio is close to but 
the magnetic ordering is collinear (with FM alignment of 
spins along the frustrated spin chains) due to the non- 
neglibible interchain couplings 1^1^ 

Based on the above discussion, we suggest that 
CuNCN undergoes long-range ordering with antiparal- 
lel alignment of spins along the c direction (due to AFM 
J) and parallel alignment of spins along the a direction 
(due to FM Ji). The ordering along the b direction 
is more difficult to find out. Both LDA and LSDA+C/ 
yield Jh < 5 K, i.e., this coupling is at the border of the 
accuracy of the present analysis. On this energy scale, 
additional factors (spin anisotropy, dipolar interactions) 
may be relevant and will influence the magnetic ordering. 
Still, it is clear that the frustrated triangular arrange- 
ment of Cu atoms within the ab plane (right panel of 
Fig. d]) does not lead to any sizable frustration, because 
the coupling is very weak. 

To get a rough estimate of the ordering temperature 
(T/v), we use theoretical results for coupled spin chains. 
Unfortunately, there are no established expressions for 
the case of different couplings along the a and b direc- 
tions. Therefore, one has to use a simplified model with 
the unique effective coupling J± between the spin chains. 
The long-range couphngs (J2, Jac) have to be neglected. 
In the following, we assume J± = (| Ji|-|- Jb)/2 ~ \ Ji\/2 — 
250 K. Similar approximations have been used for other 
spin-chain materials with spatially anisotropic interchain 
couplings and overestimated T/v by a factor of 3 — 4.^^ 
The expression by Schulz^^ igg^^^g ^ 739 k. Then, 

the actual T/v should be scaled down to 150 — 250 K due 
to the spatial anisotropy of interchain coupling. Further- 
more, one can expect an even lower T/v due to the frustra- 
tion of the spin system. Thus, the reasonable estimate is 
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100 — 150 K which is still detectable for the experimental 
methods. We suggest that the 70 K anomaly in the mag- 
netic susceptibility^i^^ should be taken as the signature 
of the long-range ordering in CuNCN. The lack of the re- 
spective specific heat anomaly and the lack of magnetic 
reflections in neutron diffraction can be explained by the 
pronounced one-dimensionality of the spin system. 

The magnitude of the magnetic neutron scattering de- 
pends on the ordered moment (sublattice magnetiza- 
tion) of the material. This quantity can be estimated 
within the two-dimensional (2D) J — Ji model, because 
2D systems have finite sublattice magnetization. Quan- 
tum Monte-Carlo simulations suggest the ordered mo- 
ment /X ~ 0.5 /xb for \ Ji\/J — 0.2.— The reduction in 
the ordered moment (compared to the classical value of 
1 /is) is caused by the strong quantum fluctuations aris- 
ing in the ID spin system. For other spin-chain materials, 
even lower /i values have been reported;-^ The reduced 
value of the ordered moment suggests that magnetic re- 
flections in neutron diffraction patterns should be weak. 
Since the proposed magnetic ordering pattern does not 
increase the unit cell (and even retains the C-centering 
symmetry in the case of FM spin alignment along the b 
direction), magnetic reflections should overlap with the 
nuclear ones, making the experimental observation of the 
magnetic ordering difficult. This may explain why the 
experimental neutron diffraction study did not indicate 
magnetic orderingfSS, while our model suggests the rather 
high Neel temperature of 100 — 150 K. To separate mag- 
netic and nuclear reflections, one can apply the polar- 
ized neutron scattering technique that has proven to be 
a sensitive tool for studying long-range ordering in frus- 
trated low-dimensional spin systems with a strongly re- 
duced sublattice magnetization.— 

The specific heat anomaly arises from the release of en- 
tropy upon the transition to the paramagnetic state. If 
the transition temperature is sufficiently small (T^r / J <C 
1), the available entropy is also small, and the anomaly 
will be completely suppressed. This effect has been stud- 
ied theoretically for square lattices with a weak interlayer 
couplings and was further confirmed by the experiments 
on Cu-based square lattice compoundsi^ For the uniform 
spin chain, a similar behavior can be expected. Assuming 
Tn/J ^ 0.030, we find S ~ 0.021i?,50 i.e., about 3 % of 
the full entropy _Rln2 for a spin-^ system or below 2 % 
of the lattice entropy at T/v (see Fig. 5 in Ref. 26) . Since 
powder samples usually show broad transition anomalies 
(see, e.g., Ref.lGll). we believe that the magnetic ordering 
in CuNCN should lead to a very weak, hardly resolvable 
specific heat anomaly at T/v- To check this hypothesis, 
specific heat measurements on a single crystal are desir- 
able. 

The above discussion shows that our spin model of 
CuNCN is consistent with all the experimental data avail- 
able so far. The 70 K susceptibility anomaly is likely an 
indication of the magnetic ordering with a low ordered 
moment that impedes the observation of magnetic scat- 
tering in conventional neutron diffraction. To study the 



magnetic ordering in CuNCN, one has to use more sen- 
sitive techniques, such as polarized neutron scattering or 
muon spin relaxation (/iSR). Resonance techniques can 
also be helpful, because they also evidence the long-range 
ordering by the shift of the absorption line in electron 
spin resonance or by the anomaly in the spin-lattice re- 
laxation rate, measured by nuclear magnetic resonance. 

The last comment regarding the experimental data 
deals with the electronic band gap Eg. LSDA-ff7 yields 
Eg = 1.0 eV for the ground-state spin configuration at 
Ud — 5 eV. This result is in qualitative agreement with 
the experimentally observed black color of CuNCN, al- 
though the resistivity measurements yield a lower activa- 
tion energy Ea ~ 0.1 eV.-^ However, one should be aware 
that resistivity measurements yield activation energy for 
the transport of charge carriers, while this energy isgen- 
erally unrelated to the electronic band gap (see Ref. |63 for 
an instructive example). In insulating transition metal 
compounds, charge is usually carried by polarons. Then, 
Ea is the lattice reorganization energy rather than the 
gap in the electronic spectrum. To obtain an experimen- 
tal estimate of Eg , optical measurements are necessary. 

Now, we will compare CuNCN to copper oxides with 
similar "edge-sharing" chains of CUX4 squares (X = N, 
O). Our results show that the introduction of the ni- 
trogen ligand preserves the physics within the structural 
[CUX2] chains. The impaired electron occupies the x^—y^ 
orbital lying within the CUX4 plaquette. Then, frus- 
trating NN and NNN interactions emerge. The NNN 
coupling J2 ~ 100 K is comparable to NNN couplings 
of 60 - 100 K in LiCusOs, LiCuV04, and other "edge- 
sharing" chain compounds. In contrast, the NN coupling 
Ji ~ —500 K is enhanced (the typical values for oxides 
are |Ji| < 200 K)^>i^>i^>^ This effect can be explained 
by the larger spatial extension of the Wannier functions 
due to the stronger hybridization between Cu and N or- 
bitals and due to the large contribution of the HOMO of 
the NCN group. Then, Hund's coupling on the nitrogen 
site leads to the strong FM interaction, as expressed by 
the Eq. (P). 

However, the above-mentioned similarities do not lead 
to the similar physics. In CuNCN, the 7r-conjugated NCN 
groups mediate the very strong superexchange coupling 
J that largely exceeds the couplings within the structural 
chains. Then, one again finds a ID magnetic behavior, 
but the spin chains run perpendicular to the structural 
chains, and the magnetic behavior resembles the "corner- 
sharing" scenario of Sr2Cu03 rather than the "edge- 
sharing" scenario of Li2Cu02, LiCu202, and LiCuV04. 
This finding reminds different alignment of structural and 
magnetic chains in vanadium oxides (e.g., Sr2V3 09^^ or 
(VO)2P207^^), although the reason is different. In vana- 
dium compounds, it is possible to suggest the correct spin 
model by a qualitative analysis based on the location of 
the magnetic orbital. In CuNCN, the in-plaquette or- 
bital favors both J, Ji, J2, and Jac (see Fig. 0]). Then, a 
careful microscopic analysis is necessary to find out the 
leading interactions and the correct spin model. 
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Despite this important difference between CuNCN and 
"edge-sharing" copper oxides, copper cyanodiimide is an 
interesting compound on its own. The very strong AFM 
couphng J ~ 2000 K along one direction is comparable 
to Sr2Cu03, only, and should lead to efficient heat trans- 
port within the spin chains. Additionally, this strong 
coupling leads to the sizable bandwidth ~ 3 eV. Such 
a bandwidth is still too small to compete with the on- 
site Coulomb repulsion. However, the application of high 
pressure can increase W and drive the system towards 
a metal-insulator transition and further unusual proper- 
ties. Experimental studies of these phenomena would be 
very interesting. 

V. CONCLUSIONS 

We have shown that copper cyanodiimide CuNCN 
should be considered as a uniform spin-^ chain system 
with spin chains running along the c direction. The inter- 
chain couplings arc found in the ab plane and show frus- 
tration within the structural [CUN2] chains. However, the 
leading interchain coupling Ji ~ —500 K is sufficiently 
strong and will likely lead to a FM spin alignment along 
the a direction. We argue that the 70 K anomaly in 
the magnetic susceptibility can be assigned to the long- 
range magnetic ordering. The ID nature of the spin sys- 



tem leads to strong quantum fluctuations that reduce the 
ordered moment down to 0.5 fis a-nd can impede the ob- 
servation of the long-range ordering in conventional neu- 
tron diffraction and in the specific heat. Further studies 
of CuNCN should include an experimental characteriza- 
tion of the magnetic ground state with polarized neutron 
scattering, /iSR, and resonance techniques. Such studies 
are currently underway or in preparation. Experimen- 
tal access to the electronic structure of CuNCN should 
be possible via optical measurements and photoemission 
spectroscopy. The potential application of this material 
could be the ballistic heat transport, while the possible 
metal-insulator transition under high pressure will be of 
fundamental interest. 
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